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Abstract 

Following the discovery of the superdeformed(SD) band in 62 Zn, we calculate several 
low- lying SD bands in 62 Zn using Relativistic Mean Field and Skyrme-Hartree-Fock 
models. Both models predict similar results, but still we can see some qualitative 
differences in the results of these two models, which are coming from the difference 
of the detail of single-particle levels. 
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Recent development of the experimental techniques enables us to study the 
nuclei at very high angular momentum. The study of the superdeformed(SD) 
bands is one of the most interesting topics in such studies. The first experi- 
mental discovery was done for 152 Dy in 1986 [1]. Since then, a large number of 
the SD bands have been observed in the A ~150, 130, 190 and 80 mass regions. 
These SD states are generated on the second minima in the potential energy 
surface which are connected with the deformed shell gaps such as N, Z ~44, 
64, 86 and 116. Based on the theoretical calculation indicating that there exist 
the N, Z ~30 shell gaps, the SD bands in the A ~60 mass region were also 
predicted[2]. In spite of much experimental effort, however, such SD states had 
not been observed due to some experimental difficulties. In 1997, a cascade of 
six 7 rays forming a new band was observed in 62 Zn[3], which was the first 
experimental discovery of the SD bands in this mass region. The extracted fa 
value is 0.45to;oy, and the SD states seem to become yrast at / ^>24 in com- 
parison with a theoretical calculation [3]. This spin value corresponds to the 
rotational frequency VL ^1.3 MeV which is the highest frequency among those 
observed so far. Besides, there are several interesting features in this mass 
region, for example, that the valence neutrons and protons occupy the same 
orbits in contrast with those in other mass regions and the residual neutron- 
proton pairing at high spin may be observed. Another example is the decay 
of these SD shapes via proton emission due to a very low Coulomb barrier. 
Such decay of the well-deformed high spin states has been already observed 
in 58 Cu[4]. For this A ~60 mass region, more and more new experimental 
data will become available in near future and so many theoretical investiga- 
tions will be done from now on, using both the semi-phenomenological models 
such as cranked Nilsson (or Woods-Saxon)-Strutinski model and fully micro- 
scopic ones, for example, cranked Hartree-Fock model with density dependent 
effective interactions and Relativistic Mean Field(RMF) model. 

RMF model is now considered as a reliable method likely to the non-relativistic 



2 



Hartree-Fock model to describe various properties of finite nuclei, not only 
/5-stable but also /3-unstable ones. Since the work of Walecka and his collabo- 
rators [5-7], RMF model has been applied to nuclear matter and the ground 
states of finite nuclei as well as the scattering data with great successes. Some 
groups also tried to apply this model to the excited states in finite nuclei. In 
1989, Munich group made a first attempt to describe the properties of rotating 
nuclei [8, 9]. They applied this model mainly to the SD bands in the A ~150, 
80 and 190 mass regions [10-13]. The RMF calculations of rotating nuclei are, 
however, limited to the nuclei in these mass regions up to now. Because RMF 
model is on its way of refinement, it is surely important to apply this model to 
various nuclear phenomena as wide as possible, and to check its applicability 
further. Therefore, in this letter, we apply RMF model to the description of 
the SD bands in 62 Zn including the newly discovered one, and examine how 
well this model can describe the very high spin states in neutron-deficient 
unstable nuclei. This is the first RMF calculation of the SD bands in this 
mass region. For comparison, we also calculate the SD bands in 62 Zn using 
Skyrme-Hartree-Fock(SHF) model. 

The starting point of RMF model is the following Lagrangian which contains 
the nucleon and several kinds of meson fields such as a-, u>- and p-mesons, 
together with the photon fields(denoted by A) mediating the Coulomb inter- 
action, 

C = Cn + At + C-uj + £p + Ca + Ant + ^NL, 

£ N =^(i 7 a d a -M)ifj, 

£ CT = i(3 a a)(dV)-im 2 a 2 , 

jC p = -^R a ^R^ + ^mlp a -p a , 

£a = -^F af3 F al3 , 
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where 

&af3 = d a UJ/3 — d/3LJ a , 

R a/ 3 = d a pp - dfsp a - g p p a X p p , 

F a /3 = daA/3 — 8/3 A a , 

are the field strength tensors. £ int is the interaction part between nucleons 
and mesons, 

Ant = - Q^T^oc 

-gpi/rfrif) ■ p a - eip~f a l ~ 2 T?l ipA a . 

In the standard applications, the non-linear self interactions among the a- 
mesons, 

r 1 3 1 4 

are also included, which are crucial for the realistic description of deformed 
nuclei. Applying the variational principle to the Lagrangian gives the equa- 
tions of motion. Within the mean field approximation, these are the Dirac 
equation for single nucleon fields ipi and the Klein-Gordon equations for the 
classical meson and photon fields. After solving these equations, we can cal- 
culate various properties of finite nuclei. 

For the application to rotating nuclei within the cranking assumption, it is 
necessary at first to write the Lagrangian in the uniformly rotating frame 
which rotates around the x-axis with a constant angular velocity Q, from which 
the equations of motion in this rotating frame can be obtained. Because the 
rotating frame is not an inertial one, a fully covariant formulation is desirable 
and we accomplished this using the technique of general relativity known as 
tetrad formalism [14, 15]. The procedure is as follows: First, according to tetrad 



4 



formalism, we can write the Lagrangian in the non-inertial frame represented 
by the metric tensor g^ u {x). Then the variational principle gives the equations 
of motion in this non-inertial frame. Finally, substituting the metric tensor of 
the uniformly rotating frame, 



/ \ 

1 -tt 2 (y 2 + z 2 ) Qz -toy 







ttz 



-Sly 



-10 



0-10 



0-1 



leads to the desired equations of motion. The resulting equations are 

| a • (-V - g u w(x)) + p(M - g a a{x)) 

+g w u°(x) - Q(L X + S^)} tpi(x) = eitpi(x), 
{-V 2 + ml-n 2 Ll}a(x)=g aPs (x), 
{-V 2 + m 2 -tfLl)u\x)=g^p v {x), 
{-V 2 + ml - n 2 (L x + S x ) 2 } u>(x)=gj v (x), 



where the p-meson and photon fields are omitted for simplicity although they 
are included in the numerical calculation. These equations are the same as 
those of Munich group. For detail, see ref.[16]. 

These equations are solved by the standard iterative diagonalization method 
using the three-dimensional harmonic oscillator eigenfunctions. Because our 
method of numerical calculation is essentially the same as that of Munich 
group [17] and its details are shown in ref. [18], we here give only the model 
parameters used. The cutoff parameters for the nucleon and the meson fields 
are taken to be A^ F =8 and A^ B =10, respectively. The parameter set called NL- 
SH is adopted which is expected to give a better description of /3-unstable 
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nuclei [19] than others. Of course the dependence on the parameter set used 
should be examined in detail, which will be discussed in a near future work. 

In the ground state of 62 Zn, all valence nucleons are in the N osc =3(pf) or- 
bits. The SD states can then be generated by putting some nucleons into the 
N osc =A intruder orbits, and these configurations are symbolically denoted as 
7i(3)~ Np (4:) Np h'(3)~ Nn (4:) Nn . According to ref.[3], the proton configurations are 
fixed to Np—2. Different SD bands are then formed depending on the number 
of neutrons lifted into the iV osc =4 orbits, which is taken as iV n =2-4 in this 
work. We consider the following SD configurations; A:7r(3)~ 2 (4) 2 z/(3)~ 2 (4) 2 , 
D:tt(3)- 2 (4) 2 z/(3)- 3 (4) 3 and C:tt(3)- 2 (4) 2 ^(3)- 4 (4) 4 . Here different configura- 
tions are possible for D according to the parity-signature quantum number 
(ir — ±,r = ±i) of the last nucleon. The occupation numbers of each parity- 
signature block for these configurations are explicitly written as 



A :tt[8 + +;8 + -; 


7- 


-+;7- 


-]i/[8 ++;8 + -; 


8- + ;8- 


-](7Ttot = + 
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7- 
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-] (TTtot = + 


r tot 


= +1), 



where the total parity and signature are also shown. The calculated dynamical 
moments of inertia of several SD bands in 62 Zn are shown in the top part of 
Fig.l. The experimental data are also shown in the figure. The sudden jumps 
in the bands A, Dl and D3 are caused by the crossing in the single neutron 
levels between the [303 7/2] (r = +i) orbit and the [312 3/2] (r = +i) one(see 
the upper part of Fig. 3). Apart from these crossings whose frequencies de- 
pend on the choice of the parameter set, the bands D3, D4 and C look to 
reproduce the experimental values well. The deformation parameters (3% are 
also calculated which are shown in the middle part of Fig.l. The experimen- 
tally extracted value corresponding to I ~24 is shown with an error bar. The 
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bands D reproduce the experimental value very well, while those of the bands 
A and C are somewhat too large or too small. From these results, we can con- 
sider that the experimentally observed one corresponds to the bands D3 or 
D4 in our calculation. The excitation energies of the bands mentioned above 
are shown in the bottom part of Fig.l relative to a rigid rotor reference. For 
/ ^25 the band A is the lowest one while the band Dl comes down for I ^25. 
The bands D3 and D4 which give the best result are located more than 1 
MeV higher than the lowest one. This result is consistent with the calculation 
of the configuration- dependent shell correction approach shown in refs.[3,20], 
but why only the band which is not the lowest one is observed experimentally 
is not clear. 

For comparison, we also perform SHF calculation of the SD bands in 62 Zn, 
using the code HFODD which was developed by Dobaczewski and Dudek[21- 
23]. The parameter set SLy4 [24] is adopted since this set is more suitable 
for /^-unstable nuclei than other sets. Besides, other sets we examined such 
as SkM* and SkP seem not to give the SD minima for 62 Zn. Figure 2 shows 
the calculated moments of inertia, the deformation parameters fa and the 
excitation energies of the same SD bands as the RMF calculation in 62 Zn. For 
the band C, the moments of inertia was already shown in Fig. 3 of ref.[25]. 
Here we can see similar trends to those of RMF shown in Fig.l, that is, the 
bands D3, D4 and C give the best results(although Dl is not so bad). As 
for f3 2 , the bands D as well as the band A can reproduce the experimental 
value. The SHF calculation with SLy4 parameter set thus also shows that the 
bands D3 and D4 are the best candidates for the experimentally observed one. 
Similarly to the RMF calculation, these bands are again located more than 
1 MeV higher than the lowest one, which can be seen in the bottom part of 
Fig. 2. So far, no parity and signature assignment for these SD bands has been 
done. If such assignment is completed and the experimentally observed band 
has negative parity, we can determine which configuration corresponds to the 
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experimental one, because the band D4 has positive signature (I n =0~, 2~, 4 - ...) 
while the band D3 has negative signature(/ 7r =l _ ,3 _ , 5 - ...). On the other hand, 
if the experimentally observed band has positive parity, both models could 
not reproduce the experimental data, which would imply that there are some 
other effects which are neglected in the present calculation, such as residual 
(neutron-proton) pairing correlations or the reflection asymmetric shapes in 
connection with the cluster-like structure. 

Looking at the results of the calculations closely, there exist some qualitative 
differences between RMF and SHF, for example, the difference of the deforma- 
tion parameters, the occurrence of the crossing, and so on. These differences 
are arising from those in the occupied single-particle orbits. Figure 3 shows 
the single neutron Routhians of the lowest SD band A as functions of the 
rotational frequency. All the orbits below the iV=32 gap are occupied. We 
notice that the position of the [310 1/2] orbit in the SHF calculation is quite 
different from that in the RMF calculation. It is located lower than the [303 
7/2] orbit, and therefore, the last two neutrons occupy the [310 1/2] orbit 
in the SHF calculation rather than the [303 7/2] orbit in RMF, which leads 
to the qualitative differences between RMF and SHF. This difference in the 
single-particle levels can be directly related to that in the deformation pa- 
rameters f3 2 which can be seen in the middle parts of Fig. 1 and Fig. 2. The 
RMF calculation gives (3 2 ~0.37(band A),0.49(band D3) and 0.46(band D4) 
at VI ~1.3 MeV, while in the SHF calculation we obtain (3 2 ~0.43(band A), 
0.50(band D3) and 0.48(band D4), respectively (both calculations give rather 
small values of 7 deformation in all of these SD bands). These values are 
consistent with the calculations shown in ref. [3], (3 2 — 0.41 — 0.49, except for 
the band A in the RMF calculation. The reason why the deformations differ 
in these two models, especially for the band A, is as follows: The [303 7/2] 
orbit is strongly upsloping with respect to /3 2 , i.e., anti-deformation-driving, 
whereas the [310 1/2] orbit is almost flat (Fig.l of ref. [26]). Because of this, 
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the RMF calculation of the band A where the [303 7/2] orbits are occupied 
gives smaller deformation than the SHF calculation in which the [303 7/2] 
orbits are empty. Both RMF and SHF calculations of the SD bands in 60 Zn, 
where the [303 7/2] orbits in the single neutron routhian are not occupied, 
give very similar deformations for the band A(/32=0.46 in both RMF and SHF 
at Q=1.3 MeV) to those of ref.[2], (3 2 =0A7 (e 2 =0.41), which strongly supports 
these discussions given above. 

The relative position of the [310 1/2] orbit and the [303 7/2] one at zero spin is 
determined by the spherical shell structure of both models. To look into this, 
we calculate the single neutron levels in the spherical core nucleus, 56 Ni. They 
are shown in Fig.4. Although the magnitudes of the L-S splittings are very 
similar in these two models, there are surely some differences, especially the 
ordering in the pf shells. In the SHF calculation, the 2py 2 and 2p 3 / 2 orbits 
are located below the I/5/2 orbit. Because of this, the N =28 shell gap reduces 
and accordingly the [310 1/2] orbit is relatively lowered in SHF in compari- 
son with the RMF calculation, which induces the different deformations and 
configurations in these two models for the SD bands in 62 Zn. 

In summary, we studied the SD bands in 62 Zn using RMF model, including the 
recently discovered one. For comparison, we also performed SHF calculation 
of the SD bands in this nucleus. In both the RMF calculation with NL-SH 
parameter set and the SHF calculation with SLy4 parameter set, the bands D3 
and D4 give the best result. There exist some qualitative differences between 
the results of these two models. They seem to be arising directly from the dif- 
ference in the equilibrium deformation, and it can be understood qualitatively 
in terms of the difference in the order of single-particle levels of the spherical 
core. Because there is only little experimental information up to now, it is 
too early to draw a definite conclusion as for whether the present models can 
describe the very high spin states in this mass region. At the present stage, 
more systematic investigation of many nuclei in this mass region using various 
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parameter sets of both RMF and SHF will be necessary, and which is under 
progress. 
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Fig.l: Moments of inertia, deformation parameters and excitation energies of 
several SD bands in 62 Zn calculated by adopting the parameter set NL-SH. 
Experimental values [3] are also shown in the top and middle figures. 

Fig. 2: Moments of inertia, deformation parameters and excitation energies of 
several SD bands in 62 Zn calculated by adopting the parameter set SLy4. 
Experimental values [3] are also shown in the top and middle figures. 

Fig. 3: Single neutron Routhians of the SD shape in 62 Zn calculated by adopting 
the parameter sets NL-SH and SLy4. The discontinuity in the upper figure is 
due to the crossing between the [303 7/2] (r = +i) orbit and the [312 3/2] (r = 
+i) one. 

Fig. 4: Single neutron levels in the ground state of the doubly magic 56 Ni 
calculated by adopting the parameter sets NL-SH and SLy4. 
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